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ABSTRACT 
This  thesis  investigates  the  problem  of  aerial  search  for  a  mobile 
surface  target.   The  initial  position  of  the  target  is  known,  and  a 
certain  amount  of  time  elapses  before  the  search  begins.   The  target 
has  a  variety  of  choices  as  to  speed  and  direction,  and  the  search  air- 
craft has  a  choice  of  search  patterns  to  be  flown.   The  method  proposed 
to  solve  the  problem  of  maximizing  or  minimizing  the  probability  of 
detection  is  the  game  theoretic  analysis  of  various  strategies  avail- 
able to  the  searcher  and  the  evader.   Four  strategies  for  each  are 
compared  and  evaluated,  with  optimal  strategies,  within  the  set  of 
strategies  investigated,  being  derived. 
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I.   Introduction. 

A  frequently  encountered  problem  of  attack  aircraft  is  the  detection 
and  destruction  of  mobile  enemy  installations  at  some  time  after  initial 
detection.   For  example,  a  flight  may  be  fired  upon  by  a  mobile  surface- 
to-air  missile  site.   Upon  return  to  the  air  base  or  carrier  the  location 
of  the  attack  is  reported  and  a  strike  against  it  is  readied  and  sent  out. 
In  the  meantime,  the  missile  site  has  often  had  time  to  depart  the  posi- 
tion in  which  it  was  originally  situated.   There  are  many  other  ways  in 
which  a  mobile  target  may  be  detected,  such  as  aerial  reconaissance,  air- 
borne electronic  detection,  or  visual  sighting.   Some  of  these  may  give 
a  probabilistic  rather  than  a  deterministic  position. 

The  mobile  target,  if  it  knows  it  has  been  discovered,  has  the  op- 
tions of  remaining  in  its  present  position  or  moving.   If  it  moves,  it 
has  some  variety  in  its  choices  of  direction  and  speed.   Often  these 
choices  of  movement  are  to  some  degree  limited.   The  search  aircraft  are 
then  confronted  with  the  problem  of  finding  and  attacking  the  moving 
strike  objective.   This  can  be  hampered  by  such  variables  as  time  to  sun- 
set, terrain,  weather,  and  other  uncontrollable  elements.   The  target  may 
be  found  by  methods  other  than  visual  if  it  is  emitting  electronic  sig- 
nals, i.e.,  if  it  has  not  moved  or  if  it  has  moved  and  is  set  up  and 
operating  in  another  position.   The  most  difficult  problem  arises  when 
the  method  of  detection  must  be  strictly  visual,  as  when  the  objective 
is  moving  or  if  it  has  set  up  in  another  position  and  is  not  emitting. 
The  degree  of  difficulty  also  depends  on  how  much  freedom  of  motion  the 
target  has.   If  the  direction  of  travel  is  limited  rather  strictly,  as 
when  there  are  only  a  small  number  of  roads  along  which  it  may  travel 
and  off-road  travel  is  not  possible,  the  searcher's  task  is  somewhat 


simplified.   If,  however,  it  is  free  to  travel  in  any  or  almost  any  direc- 
tion, the  problem  becomes  much  more  difficult.   This  thesis  examines  some 
of  the  problems  confronting  a  searcher  faced  with  the  problem  of  discover- 
ing a  target  with  unrestricted  direction  of  movement,  and  with  varying 
max  imum  s  peed  s . 

This  problem  may  be  attacked  in  two  fashions.   One  method  is  to  con- 
sider the  movements  of  the  target  as  being  determined  probabilistically, 
with  such  variables  as  direction,  velocity,  and  time  of  movement  all  being 
described  by  some  density  function,  so  that  an  estimate  of  its  most  likely 
position,  or  most  likely  set  of  positions,  can  be  found  using  the  usual 
rules  of  probability.   In  fact,  all  these  variables  have  some  probabilis- 
tic limitations,  but  within  these  limitations  there  is  another  influence: 
the  judgment  of  the  commander  of  the  site.   It  would  be  somewhat  naive  to 
think  that  the  commander  would  be  guided  solely  by  the  constraints  imposed 
by  the  state  of  nature.   Rather,  it  is  more  reasonable  to  assume  that  he 
would  try  to  vary  his  behavior  within  the  limits  of  his  physical  con- 
straints so  as  to  attain  maximum  liklihood  of  evasion.   If  this  did  not 
occur  immediately,  then  certainly  experience  would  lead  him  to  adopt  more 
sophisticated  evasive  tactics.   The  intelligent  searcher  would  then  alter 
his  search  pattern  to  counter  the  evasive  tactics  to  as  great  a  degree  as 
possible.   This  development  brings  the  problem  into  the  area  of  game 
theory,  and  it  is  with  this  aspect  that  this  paper  concerns  itself- 

To  simplify  the  problem  a  number  of  probabilistic  variables  are 
treated  deterministically.   The  effect  of  weather,  time  of  day,  terrain 
and  variance  in  terrain  over  the  search  area,  and  other  external  influ- 
ences are  ignored.   These  and  other  assumptions  are  discussed  in  the  next 
section.   A  small  number  of  strategies  for  both  searcher  and  evader  are 


discussed.   While  an  attempt  has  been  made  to  consider  the  most  advan- 
tageous strategies,  there  are  probably  many  of  merit  which  have  not  been 
included.   It  is  felt  that  an  approach  has  been  made  to  the  solution  of 
the  problem,  however,  and  some  possibly  beneficial  tactics  derived. 

The  general  approach  is  to  select  several  tactics,  both  for  the 
searcher  and  the  target.   These  tactics  are  then  evaluated  in  pairwise 
combinations,  using  game  theory  criteria.   This  allows  conclusions  to  be 
drawn  concerning  the  relative  merits  of  the  strategies. 

The  next  section  presents  the  basic  assumptions  of  the  model  and 
attempts  to  justify  them.   The  ensuing  sections  develop  the  general 
theoretical  framework  and  the  specific  mathematical  formulations  of  the 
proposed  strategies.   After  comparing  these  strategies,  some  conclusions 
as  to  the  value  of  each  are  inferred. 


II.   Theoretical  Basis  for  the  Model. 
2.1  Assumptions. 

It  would  be  an  almost  impossible  task  to  obtain  a  universal  solution 
to  this  search  problem,  in  which  all  relevant  factors  were  included.   On 
the  other  hand,  elimination  of  a  critical  variable  would  leave  the  con- 
clusions worthless.   The  assumptions  made  in  the  development  of  the  model 
are  set  forth  here  in  an  attempt  to  explain  its  limitations  and  clarify 
its  objectives. 

The  space  in  which  the  game  is  played  is  a  two-dimensional  Euclidean 
geographical  space.   The  terrain  problem  is  ignored  in  that  the  area  is 
considered  to  be  perfectly  flat,  and  in  no  way  affects  either  target  mo- 
bility or  the  searcher's  visibility.   The  target  is  free  to  move  in  any 
direction  and  at  any  velocity  up  to  its  maximum  for  any  duration  of  time. 
The  searcher  has  a  fixed  range  sighting  probability;  within  this  range  the 
probability  of  detection  is  1.0,  outside  this  range  the  probability  of 
detection  is  0.   The  searcher  is  capable  of  flying  any  pattern  exactly. 
The  initial  position  of  the  target  is  known  exactly.   The  target  knows 
he  has  been  detected,  and  takes  action  accordingly.   No  variables  such 
as  weather,  time  of  day,  or  other  external  phenomena  are  considered.   The 
time  at  which  the  aircraft  will  arrive  to  begin  his  search  is  known  to 
the  target  as  soon  as  the  initial  detection  occurs. 

There  are  other  minor  assumptions  in  the  model.   These  cited  are 
considered  to  be  the  most  relevant.   To  be  able  to  claim  that  the  model 
has  any  relevance  to  the  real  world  some  of  these  assumptions  must  be 
just  if  ied. 

First,  consider  the  playing  area.   Unless  the  terrain  is  extremely 
rough,  the  aircraft  acts  essentially  as  if  it  were  flying  above  a  two- 
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dimensional   area.      The  behavior  of  the  target    is  more   affected  by  this   as- 
sumption.     Since  the  target-chosen  factors   are   speed   and  direction,    the 
effect  of   terrain  on  these   is   relevant.      In  a  real-world   situation,    the 
target  might    indeed  be  somewhat   terrain   limited.      However,   most  missile 
and  radar   sites   are  near  cities,    which   in  turn  are  usually  situated   in 
flat   terrain.      Also,   there   is   usually  a  heavy  density  of  roads   near  any 
city,    which  gives   the  target   a  great   degree  of   flexibility   in   its  move- 
ments,  both  as   to   direction  and   speed.      If   the   area   around  the   initial 
position   is   not   travelable    in  almost   any  direction,    the  problem   is    simpli- 
fied to   some  extent,   but    is   susceptible  to  the  same   type  analysis   as   is 
proposed  here. 

The  assumption  that   the  target   can  move   at   any   speed   up  to   its  max- 
imum for  the  duration  of  the  search   is   reasonable   if  the  search  aircraft 
arrive  within  a  reasonable    length  of   time,    on  the  order  of   four  hours 
after   initial  detection,    and   stay  on   station  a  moderate   length  of  time. 
Since  most   attack  circraft   can  fly  at    low  altitudes    for  only   about   two 
hours,   the  total    time  the   target  must   be  able  to   keep  moving   is   about 
six  or  seven  hours.      This    is  well  within  the  capability  of  most  mobile 
sites. 

The  method  of   first   detection  of   the  target    is  often  visual  or  photo 
graphic,    in  which  cases   the  position   is   known  exactly.      If  the   position 
is   not   known  exactly,   the  point  of  highest   probability   is   usually  taken 
as  the  true  position  and  treated  deterministically   in  any  actual    search. 
The  effect   of    large  errors    in  the  position  of  the  original    site  is   not 
examined   in  this   discussion. 

The   implications  of   a  fixed  range  of  detection  are  well   known  and 
no   further   implications   have  arisen   in  this   analysis,     (sj   However,    the 
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idea  that  an  aircraft  can  fly  any  desired  pattern  exactly  will  strike  any 
experienced  aviator  as  a  flagrant  misconception.   In  fact,  this  is  prob- 
ably the  most  critical  assumption  of  the  model.   Recent  years  have  seen 
the  development  of  Doppler  navigational  equipment  which  enables  an  air- 
craft to  fly  a  precise  track,  and  airborne  radar  in  a  control  aircraft 
can  guide  other  aircraft  over  a  desired  pattern  with  good  accuracy.   There 
is  no  reason  to  believe  that  such  techniques  could  not  be  applied  to  search 
aircraft  if  the  situation  warranted. 

Another  critical  assumption  is  that  the  target  knows,  as  soon  as  de- 
tected, at  what  time  the  search  aircraft  will  arrive.   Some  examination 
of  the  actual  procedures  reveals  this  assumption  to  be  less  restrictive 
than  it  initially  appears.   If  the  target  moves,  the  probability  of  de- 
tection diminishes  with  time  in  any  of  the  strategies  considered.   The 
search  aircraft  will  therefore  want  to  begin  the  search  as  rapidly  as 
possible.   The  target  often  knows  the  position  of  the  base  or  carrier 
from  which  the  attack  aircraft  will  be  launched,  and  can  therefore  pre- 
dict fairly  accurately  their  arrival  time.   An  estimate  might  also  be  ob- 
tained from  observing  past  arrival  times. 

While  these  assumptions  restrict  the  application  of  the  model  to 
a  limited  situation,  they  allow  the  formulation  of  a  general  mathematical 
model.   Refinements  of  the  general  approach  may  be  made  to  deal  with 
specific  deviations. 
2.2   Outline  of  Game  Theory  Approach. 

The  searcher,  denoted  player  A,  and  target,  denoted  player  B,  have 
available  to  them  a  large  and  possibly  infinite  set  of  strategies.   Each 
attempts  to  maximize  his  payoff  by  choosing  judiciously  among  these 
strategies.   B  may  choose  any  direction  and  any  speed  up  to  his  maximum, 


12 


and  may  stop  his  movement   at   any  time.      A  may  choose  any  search  pattern 
over  the  playing  area  which  does   not  exceed  his   range.      The  payoff   func- 
tion for  A  is   the   probability  of   target   detection.      The  payoff   to   B   is 
either  the   negative  of  the   probability  of  detection,    -P(D),   or  the  prob- 
ability that  A  does   not  detect   B,    l-P(D).      If   the  former   is   defined   as 
B's   payoff   the  game   is    zero-sum;    if   the   latter  the  game   is   constant-sum 
and   the  sum   is  one. 

The  set  of   strategies    for   player  A  is   denoted  P     and  each  element  of 
this    strategy   set    is   denoted  «»<  .      The   set  of   strategies   for  player  B   is 
denoted  ^    and  the  elements   of  the  set   are  denoted  A    .      The  pairs  of 
strategies    C4*,^)    are  elements  of  the  Cartesian  product  of   Z"1   and  2L  , 
denoted  [^  x  &  .      The  payoff   to   A  for  a  given  pair  of   strategies,     "S> 
and  fl      ,    is  M(^,/J)   =   P(A  detects  B/  o<  ■  •**  1  f>  =  ($0  ) .      For  an  opti- 
mal   strategy  to   exist,   there  must   exist   strategies    c*     and    p,     such  that 

min       max     M(°<  ,  fi>  )    =  max     min     M(<*  ,  fi  )    =   M(o<   ,  A  )      [kl 
p  o<  *        ft 

whether   in  fact   such  strategies   exist  has  not   been  determined.     To 
determine  the  existance  of   such   strategies,    either   all   possible  strate- 
gies must  be  examined,   or    it  must   be   proven  topo logically  that    such  strat- 
egies  exist.      The   first  method   is  difficult   due  to   the   large   number  of 
possible   strategies.      There   exist   some  existence  theorems  which  pertain 
to   a   pursuit-and-evasion  game,    but   none  were  found  which  were  applicable 
to   the  game  under  discussion.      An  existence  theorem  developed   by  Ryll- 
Nardzewski    [l]    which  was   the  most  general   as   to   existence  criteria,    indi- 
cated that   a  solution   involving  pure   strategies   did   not   exist. 

It    is   also   possible  that   an  optimal    solution  which   is   a  convex  com- 
bination of  various   strategies   exists.      To   find   this,    it    is   necessary  to 
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examine  the  payoff  matrix  of  the  various  strategies,  and  to  apply  the 
usual  methods  of  game  theory  to  find  a  solution.   This  requires  develop- 
ing the  probability  of  detection  for  various  strategies.   Therefore  it  is 
worthwhile  to  develop  payoff  functions  for  individually  considered  strate- 
gies of  A  and  B. 

The  general  procedure  then  is  to  choose  some  strategy  oLa    and  some 
strategy  /S  .   /3  specifies  a  rule  for  choosing  direction,  speed  and 
stopping  point  for  the  target.  &(     specifies  a  pattern  to  be  flown  over 
the  playing  area.   The  payoff  is  P(A  detects  b/e<  =o<d  ***  fl  =  fi0   )•   These 
strategies  may  be  probabilistic  and  specify  a  probability  function  for 
choosing  speed  or  initial  aircraft  search  point,  or  may  be  deterministic, 
such  as  always  start  the  search  at  a  radius  of  five  miles,  or  always  re- 
main at  the  initial  position.   If  the  set  of  attractive  strategies  is 
compared  pairwise  in  this  manner,  the  matrix  of  strategy  payoffs  is  gen- 
erated. 
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III.   Development  of  Strategies. 
3.1   Strategies  Considered. 

As  has  been  noted,  there  are  a  great  many  possible  strategies  for 
both  searcher  and  target.   After  considering  a  large  number  of  strategies, 
the  following  were  selected  for  investigation  on  the  basis  of  their  intui- 
tive attractiveness  or  because  they  are  frequently  used  in  practice.   A 
glossary  of  terms  is  given  in  table  one.   The  strategies  as  proposed  by 
the  author  are: 

1 

Player  A  searches  over  v   in  the  most  efficient  manner  which  also 

searches  over  all  directions  of  possible  target  travel,  (0,  2fT").   This  is 

taken  to  be  a  spiral  search,  and  the  choice  of  the  target  velocities  to 

be  searched  is  determined  in  a  game  theoretic  manner. 

Table  I.   Glossary  of  Terms 

v   target  speed 

v   search  aircraft  speed 

Vx.    maximum  target  speed 
max 

Q      direction  of  target  movement  or  polar  angular  coordinate 

of  aircraft  position 

R   aircraft  range 

T  .   minimum  time  after  sighting  at  which  target  may  move 

T    time  from  initial  detection  of  aircraft  arrival 
arr 

t~  T    -  T  .   or  time  the  target  has  been  in  motion  at  aircraft 
0  arr    mm 

arrival  time 

r  radius  from  initial  target  position 

r  radius  at  which  the  aircraft  starts  its  search 
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S^  circumference  of  the  ith  circuit  searched  by  the  aircraft 
W   fixed  sighting  range 

A  diagram  of  a  spiral  search  is  given  in  Figure  1. 

c<   . 

2* 

This  is  essentially  the  same  as  OC  ,  with  the  additional  condi- 

3  1 

tion  that   the  original   target   position   is   always    investigated. 

3 

The  aircraft  search  randomly  in  the  circle  of  radius  (v.    *  t) , 

max 
or  the  circle  within  which  the  target  must  be. 

c*  . 

The  aircraft  search  uniformly  over  v  and  (0,  2#"  ) 
The  evasive  strategies  are 

The  target  chooses  velocity  and  direction  to  countered   accord- 
ing to  a  game  theory  criteria. 


The  target  chooses  velocity  and  direction  to  counter^   accord- 
ing to  game  theory  criteria. 

'  3' 

The  target  chooses  speed  uniformly  on  (0,  v    )  and  direction 

max 
uniformly  on  (0,  27f   ). 

The  target  chooses  direction  uniformly  on  (0,  2//~)  .      The  target 

proceeds  radially  outward  so  as  to  distribute  himself  uniformly  over  the 

geographic  area  bounded  by  r  =  v.     •  t. 

max 
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Figure   1 
Diagram  of   Spiral   Search 
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3.3  Spiral  Search. 

It  is  possible  to  demonstrate  that  a  spiral  search  for  a  target 
moving  radially  outward  from  a  fixed  point  is  the  most  efficient  way 
of  searching  over  various  target  velocities,  if  the  target  direction  is 
uniformly  distributed  on  (0,  2/7*).      This  may  be  intuitively  justified  as 
follows:   Suppose  a  searcher  has  a  given  fixed  range.   He  desires  a 
search  pattern  which  will  minimize  his  path  length  while  not  searching 
over  any  velocity  range  more  than  once.   For  a  single  circuit,  i.e.,  a 
search  of  9  from  0  to  2^,  the  target  velocities  depend  upon  search  width 
W  and  the  initial  radius  of  the  search,  rQ.   If  the  velocity  range  is  to 
remain  constant  throughout  the  circuit,  the  search  radius  must  move  out- 
ward with  the  velocity  which  would  have  put  the  target  at  r  ,  which  is 
r^/t^.   If  the  searcher  moves  outward  in  such  a  manner,  and  increases  the 
radial  component  of  his  velocity  with  each  subsequent  circuit  such  that 

r  «  r   n  -  (r   .  /  t~)  •  t,  he  will  not  re-search  any  previously  searched 
n    n-l     n-1    U 

velocity.   This  pattern  may  be  flown  in  several  ways.   The  outward  radial 
component  of  velocity  may  be  chosen  to  remain  uniform  throughout  the 
search.   Another  method  is  to  vary  the  radial  component  after  the  comple- 
tion of  each  circuit,  thereby  flying  a  series  of  spirals  joined  by  short 
cross-legs.  While  this  is  slightly  less  efficient  than  a  continuous 
spiral,  it  has  the  advantage  that  it  searches  over  a  continuous  target 
velocity  range.   This  makes  it  more  compatible  with  an  analysis  of  the 
optimum  velocity  choices  of  the  target,  and  for  this  reason  this  method 
is  used. 

The  equation  for  a  single  circuit  is  derived  as  follows: 
The  search  aircraft  flies  an  increment  of  the  circuit  of  length  A  S 
in  time  A  t.   During  this  time  the  target  moves  radially  outward  with 
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r  =  v.»<*  t  =  v 


velocity  v. 


or  for  a  complete  circuit,    r  = 

To  derive  S,  take  9  to  be  in  radians. 
(1)     dS  =  rd6 


vt 
then  since  r  =  — _  S, 

v 


**< 

since 

t   = 

S 

V 

a 

^      s 

va 

a 

v 

v 


— —  sde 


or 


a 


H_  =  J[t_  de 


va 


Integrating, 


i+-R 


vt 
d9  ^       In  S  +  cn  =  — —  0  +  Co 

1     v_         2 


or  equivalently, 

(2)     S  =  Cj  e  a   +  c2 


Since  (1)  r  = 


d9~  ' 

Vt  9 

(3)     r  =  —1-  c,e  a 
v_    1 


Imposing  boundary  conditions, 

r  =  rc 
6  =  0 

S  =  0 


va 

c,  =  r   from  (3) 

1    vt   o 


v 

c0  = r   from  (2) 

z  v     o 

t 
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The  final  forms  of  the  expressions  are 


S   =  — *-  r      (   e  Va        -   1    ) 
v+.        o 


vt 

r   =   rQ   e     a 


v 


t 


—  e 

V 

To  get   these  expressions    in  a  more  manageable  form,    express   e      a 

as   an   infinite  series. 

2 
vt  v.  o 

e  Va        =   1    +     Vt     9   +     v^  +    .    .    .    . 

— 
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For  v    »  v    , 

r  a 


lie 

e  va     ^     1    ♦  _£_  9 

V 

a 


Since  v     is   usually  over    300  and  v      is   usually   less   than   30,    this 

a.  L 

is  a  valid  approximation. 


Then 


va  vt 

S  =  — —  r   (  1  +  — —  9  -  1  )  =  9r 
vt    o       va  o 


and 


vt 
r  =  r   (  1  +  — —  9  ) 

° 
At  the  end  of  each  circuit,  the  searcher  must  position  himself  to 
search  over  a  new  velocity  increment.   If  the  detection  range  either  side 
of  the  flight  path  is  W,  the  minimum  velocity  unsearched  is  (r~  +  W)/  t  , 
which  is  that  velocity  which  would  put  the  target  at  r.  +  W  at  the  start 
of  the  search.   Since  the  search  width  is  2W,  a  radius  of  rn  +  2W  will 
search  velocities  from  (tq  +  W)  /  t0  to  (r0  +  3W)  /  tg,  if  the  velocity 
is  assumed  to  be  constant  across  the  search  width.   This  is  not  strictly 
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F0  +  W   .   r0  +  2w   ^   r0  +  3W    u        •     o.u-   j- 
true,  since  ^  <£  _  ;  however  since  this  dis- 

t0         to         t0 
tance  is  small  it  will  be  ignored.   For  example,  if  W  =  1,  t„  =  2,  then 
W/t^  is  .5,  which  is  the  maximum  deviation  from  the  velocity  at  the  cen- 
ter of  the  search  strip.   As  the  time  increases,  this  deviation  becomes 
smaller. 

If  the  target  were  at  Tq  +  2W  at  tQ,  it  would  move  radially  outward 
during  the  time  the  first  circuit  was  being  flown.  The  time  required  to 
search  the  first  velocity  increment  is 

S0       2   r0 


v„        v 

a         a 

Then  the   initial   radius   for  the  following  circuit    is 

rQ    +   2w  SQ 

rn    =   rn   +   2W  +  -2 (—1-) 

10  t0 

The  distance   from  the  position  at   the  end  of  the   first   circuit   to   the 
beginning   point  of  the  next   circuit    is 

2     r0                            4W     rn 
r   =   r,    -  rn    (    1   +  —  )    =   2W  +  iL. 

1  °  Vo  ^a 

If  r     =   100,    W  =    1,   va   =    300,    and  tQ   =  2, 

r   =   2   +      8      -10°      =   12.16 
2-300 

Since  these  cross-legs   are   small   compared  to   total    aircraft   range,    they 

will   be   ignored    in  further   computations. 

Then  to   compute  any   r,    it    is   necessary  to   take   its   position  at   t0 

and  to   add  to    it  the  distance  traveled  while  the   aircraft   was    searching 

n-1 

previous    circuits.  ^~ 

<_  s. 

rn   +   2nW  4-ci      x 

r=   rn   +   2nW  +    (— - )      -i-± 

u  t~  v 

^0  a 
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since 

S  =  2^r  , 


n-1 


S  =  2 
n 


T  rn  +  2nW     .  n  Si   / 


n-1 
Because  of  the  term  ^_  S.,  a  recursive  relation  exists. 

i=0      X 

A  solution  for  N   =   ri     „   is   possible   if   the  total   path  flown  is   re- 

max    r  r 

stricted  by  the  aircraft  range,  or 

N 


R  =  2-  SL 

i=0 
However,  because  of  the  complexity  of  the  expression  an  analytic  solution 

is  very  difficult  and  an  iterative  method  of  solution  for  the  maximum 

velocity  searched  as  a  function  of  Tq  is  laborious.   For  these  reasons  a 

computer  program  was  used  to  solve  the  expression.   Thus  a  method  exists 

for  solving  for  the  range  of  target  velocities  searched  for  any  r  ,  with 

input  parameters  W,  R,  v  ,  v.     and  tn.   This  solution  is  necessarily 

a        rmax  u 

discrete,    since  the  Tq's   are  discrete;    however,    any  degree  of  accuracy 

desired    is   obtainable. 

3.4     Game  Theory  Evaluations  of  the  Proposed  Strategies. 

3.4.1     Discussion  of  StrategyoC      . 

Using  the   previously  derived  expression  for  the   range  of  velocities 
searched  over   for  any   r    ,    the  minimum  and  maximum  target  velocities  may 
be  obtained.      The  probability  of  detection  for  any  given   initial   search 


radius,    say   r    ,    is 


P  f  v    .      (r*)  <    v    <    v  (r*)    /  r      =  r*  / 

t    mm       o  t  max       cr    /      o  o/ 
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Then 

M(°V   '31)    =  P  ^ 

=  P  (v   .      (r*)  <   v    <  v  (r*)  I  r     =  r*Jp/r     =  r*  J 

(    mm       o  t  max       cr  /     0  o/       {   o  o  f 

The   lower  the  value  of  r     the  greater   is   the   range  of   target   vel- 
ocities   searched.      For  the   last   circuit,    generally  not    all   values   of   9 
are  searched   since  the  circuit  must  terminate  when  maximum  range    is   reached. 
The  target   direction  of  travel    is  distributed  uniformaly  on   (0,    2      ),   which 
means   the   probability  of  detection  on  the   last   circuit    is  the   percentage 
of  the  circumference  searched.      This  probability   is   denoted  9N-      The  den- 
sity  function  for  P(D/rQ   =   r~)    is 


r0-w 


0  v    /    =  v    .       (r   ) 

t^  t  min   ^r(T 


rQ    -  W       .  rQ    +   2NW-   W 

-    v.  < 


f<#    =        <v  *0  '  fc0 


rQ    +    2NW  -   W  rQ    +    2NW  +  W 

<    vt<    =   v        (r   ) 

t  ^      c>*  t  max     0 

r0  0 


3  v  (r*) ^  v 

max        0    ^      t 

The  graph  of   this  density   function   is   depicted   in  Figure   2. 

Denote  the  range  of  velocities   searched  by  v   .      The  payoff   function 

is  M    (v    .    v)=P(D/v     =v      and  v.     =   v.) . 
r7      t  r  r  t  t 


The  distribution  functions  for  v  and  v.  are  F(v  )  and  G(v.  ).   These 

r  t  r  t 

are  not   fixed  distributions,    but   are  any  possible  distributions.      The  ex- 
pected value  of  the   game   for  two  distributions   F  and  G   is   denoted  E(F,G). 

vt  vt 

rmax      .-'•max 
M(vr,    vfc)    dF(vr)    dG(vt) 

0  0 
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Pfl>) 


i.ft 


H, 


rQ4aA/»*-»/ 


r^-f  a/UM»w< 


^r 


Figure  2 
Conditional  Density  Function  of  the  Probability  of  Detection 

If  M(v  ,  v  )  is  continuous  then  a  solution  is  guaranteed.   However,  in 
the  case  of  spiral  search  M  is  not  continuous,  and  therefore  no  solution 
is  guaranteed.   However,  an  examination  of  the  proof  of  this  theorem  in 

4   indicates  that  the  requirement  of  a  continuous  payoff  function  may  be 
too  restrictive.   The  assumption  of  continuity  is  necessary  in  the  proof 
to  guarantee  the  existance  of  the  Stieltjes  integral  of 

J  M(x,  y)  dF(x) 


and 


; 


M(x,  y)  dG(y) 
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and  to  guarantee  the  existance  of  certain  minima  and  maxima  over  the  do- 
main of  the  payoff  function.   It  has  not  been  determined  if  the  existence 
of  a  solution  to  the  problem  under  discussion  exists. 

The  existence  of  a  game  theory  solution  has  been  presumed  and  den- 
sity functions  for  the  searcher  and  target  derived.   Due  to  the  complex- 
ity of  the  analytic  expression  giving  the  range  of  target  velocities 
searched  for  a  given  rn,  an  approximation  to  the  solution  using  discrete 
values  of  rn  and  v.  has  been  used. 

If  rn  and  v  are  allowed  to  take  on  discrete  values  and  the  probabil- 
ity of  detection  derived  for  each  combination  of  these  discrete  values,  a 
matrix  approximating  the  continuous-domained  payoff  function  may  be  ob- 
tained.  The  degree  of  accuracy  is  dependent  on  the  fineness  of  the  incre- 
ments of  v  and  r  ;  for  purposes  of  the  present  investigation  increments 
of  one  unit  of  speed,  i.e.,  one  knot  or  one  mile  per  hour  were  used.   The 

resultant  matrix  is  m  x  m+1,  where  m  =  v.     .   There  are  several  meth- 

max 

ods  to  obtain  the  game  theory  solution  of  this  matrix,  but  the  precise 
methods  of  solution  have  the  detriment  of  being  extremely  complex  for  a 
matrix  of  any  size.   For  this  reason  an  algorithm  which  gives  an  approx- 
imate solution  was  used.   The  algorithm  and  the  computer  program  used  to 
perform  the  calculations  are  discussed  in  Appendix  I.   An  extensive  dis- 
cussion with  proof  of  convergence  is  given  in  Karlin  f2j  .   The  relevant 
properties  of  the  method  are: 

(1)  An  approximate  strategy  for  A  and  B  may  be  found. 

(2)  Upper  and  lower  bounds  on  the  value  of  the  game  are  generated. 
The  accuracy  of  the  solution  depends  on  the  number  of  iterations  per- 
formed.  Two  thousand  iterations  are  considered  to  give  a  usable  approx- 
imation for  the  matrix  under  discussion.   As  an  accuracy  check  one  program 
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was  run  which  performed  twenty  thousand  iterations.   The  upper  and  lower 
values  of  the  game,  denoted  V  and  V,  were  .1254  and  .0845  for  2000  iter- 
ations and  .1153  and  .1017  for  20,000  iterations.   The  total  computer 
time  time  for  the  more  accurate  method  is  four  minutes  and  seven  seconds. 
With  2,000  iterations  ten  sets  of  input  parameters  may  be  run  in  three 
minutes  and  fifty  seconds.   Using  this  program  optimal  strategies  for  A 
and  B  were  found,  with  A  restricted  to  using  only  spiral  search.   The 

strategies  for  searcher  and  target  for  varying  v     and  t   are  tabulated 

max      ° 

in  Appendix  II. 

3.4.2.   Discussion  of  Strategy  o(     . 

All  of  the  theory  and  general  discussion  relating  to  oC   is  appli- 
cable too<  J*   The  same  computer  program  with  appropriate  modifications 
was  used  to  solve  the  payoff  matrix.   The  range  of  the  search  aircraft 
was  diminished  by  the  radius  of  initial  search  to  compensate  for  search- 
ing v  =  0  on  each  search  pattern,  which  had  little  effect  on  the  out- 
come.  This  strategy  forced  B  to  always  move  from  the  position  of  initial 
detection,  which  to  some  extent  alleviates  the  assumption  that  the  time 
the  target  begins  its  movement  is  a  deterministic  rather  than  a  probabil- 
istic quantity.   The  reason  for  this  assumption  is  that  there  is  a  lower 
bound  on  the  length  of  time  the  site  takes  to  disassemble  components  and 
makes  ready  to  move.   If  B  is  forced  to  move,  the  longer  he  is  in  motion 
the  more  difficult  is  the  problem  of  detection.   If  B  knows  the  initial 
position  will  always  be  searched,  he  is  more  likely  to  make  an  effort  to 
become  highly  efficient  in  preparing  to  move,  and  therefore  the  range  of 
times  at  which  he  moves  cluster  more  closely  around  the  minimum  move  time. 
The  solution  generated  by  the  computer  for  this  strategy  indicates  that  B 
will  move  at  a  very  low  speed  in  about  the  same  proportion  that  he  did  not 
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move  when  A  plays  strategy  ^ , . 

3.4.3.  Discussion  of  Strategy©^  . 

This  strategy  requires  A  to  search  randomly  in  a  circle  of  radius 

(tQ  "v    ) .   It  has  the  advantage  of  not  requiring  any  expenditure  of 
max 

effort  on  navigation,  and  in  fact  is  often  used  in  practice.  It  is  an 
inefficient  search  method  and  is  included  largely  as  a  basis  of  compari- 
son for  more  effective  strategies.  The  area  in  which  it  searches  is  not 
constant  if  the  target  is  moving  radially  outward,  and  the  probability 
of  detection  is  less  than  is  given  in  the  following  formula.  However, 
since  the  probability  of  detection  is  low  for  this  method,  no  accuracy 
was  lost  in  the  general  conclusions. 

The  probability  of  detection  is  given  by 

-2WR 
PCD)  =  1  -  e   A         £} 

The  derivation  of  this   formula   is  given  by  Koopman.     f~3J 

3.4.4.  Discussion  of  Strategy0*    . 

In  this  strategy  A  searches  all  velocities  of  B  with  equal  likeli- 
hood, or  v  is  distributed  uniformaly  on  (0,  v    ).   There  is  a  problem 

max 
in  this  strategy  arising  from  the  fact  that  searching  some  velocities 

requires  less  distance  traveled  than  others.   The  strategy  is  therefore 
formulated  as  the  sequential  choice  of  velocities  to  be  searched  with 
each  choice  being  uniformly  distributed  on  the  target  velocity  range. 
The  expected  probability  of  detection  based  on  the  expected  range  of 
velocities  searched  is  used  in  the  payoff  matrix  of  the  various  strate- 
gies. 

3.4.5.  Discussion  of  Strategies  /5   and  A     . 

These  strategies  are  defined  in  the  course  of  deriving  strategies 
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°*.  and  °^o-   As  with  ^ -,  and  O^    discrete  values  are  used  to  approximate 
the  continuous  strategy. 

3.4.6.  Discussion  of  Strategy  /3  . 

B  moves  at  time  t~  with  a  velocity  chosen  uniformly  on  (0,  v.    ) 

max 

and  continues  to  move  with  this  velocity  throughout  the  time  period.   The 

direction  of  travel  is  chosen  uniformly  on  (0,  2lt  ) .   This  distribution 

is  also  obtained  if  the  target  moves  with  velocity  v     and  stops  at 

max 
times  uniformly  distributed  on  (0,  t^)  where  0  denotes  the  minimum  time 

at  which  the  target  may  move. 

3.4.7.  Discussion  of   Strategy/3 

In  this   strategy  player  B  distributes   himself   uniformly  over  the   geo- 
graphical   area  of  radius  v  .    tn.      The  probability  of  detection  is 


max 


where 


PCD)    =  -£- 
A 


a  =  7Tr2 


The  distribution   function   in  terms  of  radius,   r,    is 

2  2  2 

F(r)    =   P(D)    =  r  =  -. r     -    .u        = T 

A  ™t V  <vt  V 

The  density  function  is 

f(r)  =  dF(r)  =  — _-.  dr 

(vt    tQ)2 
max  u 


Since  r  =  tgV, 


dr  =  tndv 


The  density  function  in  terms  of  v  is 


2tnvt~dv       0 
cr   \               °   Q         2v     , 
f(v)  =  n 7—   =  tj-      dv 


m 


m 
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and   the  distribution   is 


P(vt  <   vr)    = 


V  9 

f  (v)   =  2- 

'  0 


Since  strategies  O^       and  oi      have  been  approximated  using  discrete  vel- 
ocity increments  this  strategy  may  be  defined  in  velocity  increments  using 
the  distribution  function,  which  gives  an  expression  compatible  with** 
ando<   . 
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IV.   Evaluation  of  Strategies. 

4.1.   Method  of  Evaluation  of  Strategies 

The  pairwise  comparison  of  all  combinations  of  strategies  considered 
may  be  made  and  the  results  , tabulated  in  a  matrix  of  payoffs.   Examina- 
tion of  this  matrix  then  reveals  the  optimum  or  combination  of  optimum 
strategies.   The  elements  of  the  matrix  are 


i      =p[D|*  =  o<in    ft  =/3    ] 


J* 

The  general   expression  for  the  probability  of  detection  for   strategies 


at    ..    and  o<  2    is 


vt  vt 

max  max 


PCD)    =     £  P(rQ   =   r.)    P(D  /  rQ   =  r .)    =  £_   "   PCrQ   =   r.)^ 

i=l  i=l 

A>Cv    .      Cr.)        v.        v  Cr.)) 

V.        mm        x  t  max        l 

The  comparison  of  o(  ,    and  (Z     ,    and  of  o*  ~   and  l3       have   been  made   in  the 

solution  of  the   game.      The   solution  for  elements   a   _    and   a91    may  be  made 

using  the  general    formula  directly.      To   solve  for  elements   a      ,    a      ,    a      , 

a24  the  following  modification  of  the  basic   formula   is   used: 


v<.  v  Cr.) 

tm„„  max        is 


'-max  "^A        *• 

PCD)    =   £2        p(r0   =   ri)       J  f(vt)  for  fi  3  and  fi 


4 


where 


i=1  v    .      Cr.) 

mm        x 


v  9  2 

max  vz        -  v 


; 


~,      s            max          mm  e        a 

f  Cvt)    =  2 for  r 


v. 


t 


4 


mm  max 


and 


v 
max 


v  -  v 

max  mm 


;max           mm  c        /» 

f(vt)=— ~ for/3 


v    •  max 

mm 
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For  0(       the  formula  for  the  probability  of  detection  of  a  uniformly  dis- 
tributed target  is  used  although  this  is  not  true  for  strategies  p  , 
(*  ~  and  £>     .   However,  the  low  probability  of  detection  for  random  search 
indicates  that  more  precise  computation  is  not  necessary. 

For  /S  ^  there  is  a  difficulty  in  that  path  length  for  searching  dif- 
ferent velocities  is  not  constant.   This  has  been  handled  by  using  the 
expected  value  of  the  path  length  for  a  single  circuit.   The  expected 
value  of  v  is  calculated  for  each  strategy  of  B  and  this  is  used  to  get 
the  expected  number  of  circuits.   This  is  then  used  to  calculate  the  prob- 
ability of  detection. 

The  expected  target  velocity  is 

vt 
max 

E(v  )  =£.    v    P(v.  =  v   ) 
t  ui     r    ri 

0 


Then 


Thus 


E  fpath  length]   =  e[L]   =2  fttQ   E(v  ) 


E  ^number  of  circuits}   =  EfcJ  =  — ■  - 

While  it  would  be  possible  to  solve  for  the  exact  probability  of  detec- 
tion of  oi  .     for  (1       and  />  y   this  was  not  done  since  rough  calculations 
indicated  that  such  a  value  would  be  very  low.   The  expected  number  of 
circuits  may  be  used  to  solve  for  a^o  directly  by 

2W 
E(C)  x  tQ 


vt 


max 


since  the  target  speed  is  distributed  uniformaly  on  (0,  v.    ).   The  ex- 


max 


pected  probability  of  detection  for  a^  has  been  calculated  to  be 


30 


2/2  tQ 

vt 
max 

for  a  single  circuit.   Therefore  a.,  was  taken  to  be 

mm      W 

2  yr  t7 

— ! _   (e(o) 


vt 

max 

The  payoff  matrix  for  parameters  R  =  500,  tQ  =  2 ,  W  ■  1 ,  v    =30 

max 
and  v2  =  300  is 


1 

2 

3 

4 

.108 

.103 

.180 

.071 

.103 

.101 

.174 

.074 

.041 

.041 

.041 

.041 

.068 

.064 

.077 

.054 

1 

2 
3 

4 

This  matrix  has  a  saddle  point   at   a„^.     Therefore  pure  strategies  exist 
for  players  A  and  B.      It   appears  that   searching  by  strategies  ©C  .    and 
0^2   i-s  highly  desirable  for  player  A,    since  the  values  of  a^  .  and  a?,- 
are  very  close  and  a  deviation  by  B  from  strategy  f}>       greatly  increases 
the  probability  of  detection.      If  B  distributes  himself  uniformly  over 
his  range  of  velocities  then  A  increases  the  probability  of  detection 
by  more  than  a  factor  of  two   if   strategies  «<  ^   and  0<  2   are  used.      To 
determine  whether  the  existence  of  a  saddle  point   in  this   case  was   unique, 
values  of  the  elements  of  the  first  two  rows  of  the  matrix  were  calcu- 
lated for  six  other   sets  of   input   parameters   and   in  each  case  a  saddle 
point   existed  at  either  a^  or  a... 
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V.      Conclusions. 

5.11      Evaluations  of   Strategies 

The   low  value  of   strategies        „  and 
vfc   searches,    indicate  that  their   use   in  an  actual   situation  would  not   be 
optimal.      However,    the  use  of   a   simple   strategy  that   corrects   some  of  the 
deficiencies  of  these  strategies  may  be  used.      If   an  area  search   is   used 
which  does   not   re-search  any  area,    the  probability  of  detection  will   be 
approximately 

2WR 

vt        *    fc0 
max 

This  is  very  slightly  less  than  the  probability  of  detection  using  strat- 
egies °^-  and  °*p,  the  spiral  searches,  and  is  very  much  simpler  to  apply. 

It  must  be  concluded  that  a  uniform  geographic  distribution  of  the 
target  is  optimal  in  the  set  of  strategies  discussed.   Since  no  counter 
to  this  strategy  seems  to  exist,  this  may  well  be  an  optimal  solution 
over  all  strategies  available  to  player  B.   If  the  equation  for  probabil- 
ity of  detection  used  for  non-overlapping  flight  path  over  a  circle  of 

radius  v    •  t   is  used,  a  very  close  approximation  to  the  probability 

max 
of  detection  given  by  the  saddle  point  is  obtained.   This  indicates  that 

strategies  c*  -■  and  o<  n   are  optimal  in  the  set  of  strategies  studied  by 
virtue  of  their  characteristic  of  non-overlapping  flight  paths  rather 
than  because  they  derive  an  advantage  because  of  some  other  characteris- 
tic.  If  this  is  true  and  if  B  plays  strategyl3  ^,  uniform  geographical 
distribution,  the  effort  needed  to  apply  strategies^  ,  and  ©<   seems 
unjustified.   However,  the  distinction  between  uniform  geographical  dis- 
tribution and  uniform  velocity  distribution  is  not  obvious,  and  conceiv- 
ably a  naive  player  B  could  choose  to  play  strategy  ($>       and  distribute 
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uniformly  on  (0,  tt    )•   Since  using  strategies  &t      or  <*(      in  this  case 

max 

would   considerably   increase  the  probability  of   detecting  the  target,    the 
anticipated   benefits  might  outweigh  the   increased  effort.      This  might   be 
augmented  by  an  attempt   to  estimate  the  type  of   strategy  applied  by  the 
evader  by  goodness-of-f it   tests,    to  determine  the  probability  that  B   is 
using  some   strategy  other  than  ^  ^  and  to   use  a  counter-strategy  which 
takes  maximum  advantage  of  B's  mistakes. 
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APPENDIX  I 
Computer  Program  for  Solving  for  Strategies    ..  and   2 
This  program  has  two  basic  parts.   The  first  develops  the  range  of 
target  velocities  searched  for  various  values  of  rn  using  a  spiral  search, 
and  generates  a  payoff  matrix  using  increments  of  one  unit  of  speed  be- 
tween successive  values  of  Tq  and  vfc.   The  second  part  of  the  program 
solves  the  matrix  using  a  method  for  approximating  the  solution  of  matrix 
games . 

The  inputs  are: 

RO       initial  radius  searched 
VA       search  aircraft  speed 

TO       length  of  time  the  target  has  been  in  motion  when 
search  begins 

RNG      search  aircraft  range 

DELR     amount  the  initial  radius  is  incremented  for  each 
successive  search 

W        search  width 

TMAX     maximum  target  speed 

The  outputs  are: 

VMIN     minimum  target  speed  searched  for  a  given  RO 

VMAX     maximum  target  speed  searched  for  a  given  RO 

THETA    percent  of  circle  searched  on  last  circuit 

RO       radius  corresponding  to  VMIN,  VMAX  and  THETA  on  same 
line  of  printout 

A        matrix  of  detection  probabilities 

X(i)     approximation  of  optimum  percent  of  time  the  search 
should  start  at  RO  equal  to  i  •  t«  +  W 

Y(i)     optimum  percent  of  time  the  target  should  choose 
vt  =  i  -  1 
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VI  upper  value  of  game 

V2  lower  value  of  game 
Internal  symbols  are: 

SUM  total  distance  flown 

S  circumference  of  present  circuit 

SLAST  percent  of  circumference  searched  on  last  circuit 

VT(j)  target  velocity  equal  to  (j-1) 

vt 
umax 

SUMA  (IA)   j£    a^  •  for  j    .}  SUMB(j)  is  maximum  over  j 
i=l 


vt 
umax 

+  1 

SUMB(IB) 

j  =  l 

a.  . 

NA 

minimum 

SUMB  (j) 

NB 

max  imum 

SUMA  (i) 

a.  .  for  i   . ^  SUMAC i)  is  minimum  over  i 
1  >  J 


YB  (IE)   number  of  times  the  target  has  picked  velocity  YB 
(IE  -  1) 

XA  (IX)   number  of  times  the  search  aircraft  has  picked  initial 
velocity  searched  XA—  IX  up  to  present  iteration. 

This  program  is  written  in  Fortran  63  for  use  on  the  CDC  1604  com- 
puter.  The  first  part  of  the  program  uses  the  equation  developed  in 
Section  3.3  to  generate  minimum  and  maximum  target  velocities  searched 
for  various  values  of  RO .   The  initial  RO  is  an  input  to  the  program,  as 
is  the  increment  of  RO,  designated  DELR.   For  purposes  of  this  investi- 
gation increments  of  one  unit  of  speed  were  desired,  and  therefore  the 
statement  DELR  =  TO  has  been  inserted.   RO  normally  begins  at  W,  which 
assures  search  of  the  initial  target  position.   Each  successive  RO  de- 
fines a  row  of  the  payoff  matrix.   After  VMIN  and  VMAX  are  found,  the 
spectrum  of  target  velocities  are  subjected  to  IF  tests  to  find  if  they 
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fall  within  or  without  the  detected  velocities.   They  are  assigned  prob- 
abilities of  detection  of  0  or  1.   Then  the  element  of  the  row  correspon- 
ding to  VMAX  is  assigned  a  probability  of  detection  THETA.   This  process 
is  continued  until  RO  =  VTMAX  ■  TO.   At  this  point  the  matrix  is  com- 
plete. 

Part  II  of  the  program  obtains  an  approximate  game  theory  solution 
for  t  is  matrix.   Obtaining  an  exact  solution  for  a  large  matrix  is  diffi- 
cult unless  a  saddle  point  exists,  which  requires  that 

max  min  a-  -  =  min  max  a-  . 

But  for  each  i  in  the  matrix  there  exists  an  element  a.  .  such  that 
ajj  =  0.   For  each  i  there  exists  an  a.  .  such  that  a.  .   0.   Therefore 
no  saddle  point  exists. 

The  exact  methods  of  computation  involve  examination  of  all  adjoints 
of  the  submatrices  of  the  matrix.   The  generation  of  the  adjoints  re- 
quires evaluation  of  the  inverse  and  the  determinant  of  the  submatrix. 

30, 
The  number  of  submatrices  of  a  30  x  30  matrix  is  Z~      (^  ) (^  ) .   There- 

n=l 

fore  a  method  approximating  the  actual  solution  has  been  used.   This  uses 
an  iterative  process  in  which  player  A  chooses  his  strategy  by  looking 
at  the  column  vector  which  is  the  sum  of  all  columns  chosen  by  B,  and 
playing  the  row  corresponding  to  the  maximum  element  in  this  vector.   B 
looks  at  a  row  vector  composed  of  the  sum  of  all  rows  chosen  by  A  and 
chooses  the  minimum  element  for  his  next  play.   The  number  of  times  each 
strategy  is  chosen  divided  by  the  total  number  is  an  approximation  of 
the  optimal  strategies.   The  minimum  element  in  the  row  vector  and  the 
maximum  element  in  the  column  vector  divided  by  the  number  of  iterations 
give  lower  and  upper  bounds  on  the  value  of  the  game. 
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The  row  vector  of  sums  is  SUMB  in  the  program,  and  the  column  vector 
of  sums  is  SUMA.   The  number  of  times  each  row  and  column  has  been  chosen 
are  XA  and  YB.   Statements  70  through  80  choose  the  maximum  over  SUMA. 
Statements  74  through  77  choose  the  minimum  over  SUMB. 

The  program  as  presented  generates  the  solution  for  strategy^  . 
To  obtain  the  matrix  and  solution  for  strategy °i   ,  ,  statement  97  and  the 
preceding  statement  may  be  deleted,  along  with  the  fourth  statement  after 
statement  17. 
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